An Instance of Failure for the MATLAB Explicit ODE45 Solver
نویسنده
چکیده
We consider the adaptive strategies applicable to a simple model describing the phase lock of two coupled oscillators. This model has been used to show an instance of failure of the ODE45 RungeKutta-Felberg solver implemented within the MATLAB ODE suite, see [J. D. Skufca. Analysis still matters: a surprising instance of failure of Runge-KuttaFelberg ODE solvers. SIAM Review, 46:729-737, 2004]. We compare the numerical results obtained with: the MATLAB ODE suite’s explicit solvers, and the local linearity strategy implemented with the classical fourth-order Runge-Kutta method as a basic method.
منابع مشابه
Steady-State and Dynamic Simulations of Gas Absorption Column Using MATLAB and SIMULINK
Separation is one of the most important process in all the chemical industries and the gas absorption is the simplest example of separation process which is generally used for the absorption of dilute components from a gaseous mixture. In the present work, a dynamic system of mathematical equation (algebraic and differential) is modeled to predict the behavior of the absorption column using mat...
متن کاملVectorized Solution of ODEs in Matlab
Vectorization is very important to the efficiency of computation in the popular problem-solving environment Matlab. It is shown that a class of Runge-Kutta methods investigated by Milne and Rosser that compute a block of new values at each step are well-suited to vectorization. Local error estimates and continuous extensions that require no additional function evaluations are derived. A (7,8) p...
متن کاملExplicit parallel two-step peer methods
The construction of s-stage explicit parallel two-step peer methods for the solution of non-stiff initial value problems with optimal zero stability and order p = s is considered. By an additional condition superconvergence of order p = s + 1 can be achieved. Further criteria for good methods are large stability regions and small error constants. Corresponding parameter sets have been obtained ...
متن کاملA comparison of solution methods for Mandelbrot-like polynomials
We compare two different root-finding methods, eigenvalue methods and homotopy methods, using three test problems: Mandelbrot polynomials, Fibonacci-Mandelbrot polynomials, and Narayana-Mandelbrot polynomials. For the eigenvalue methods, using both MATLAB and Maple, we computed the eigenvalues of a specialized recursively-constructed, supersparse, upper Hessenberg matrix, inspired by Piers Lawr...
متن کاملCoupling Nonlinear Element Free Galerkin and Linear Galerkin Finite Volume Solver for 2D Modeling of Local Plasticity in Structural Material
This paper introduces a computational strategy to collaboratively develop the Galerkin Finite Volume Method (GFVM) as one of the most straightforward and efficient explicit numerical methods to solve structural problems encountering material nonlinearity in a small limited area, while the remainder of the domain represents a linear elastic behavior. In this regard, the Element Free Galerkin met...
متن کامل